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Abstract 

Cauchy Biorthogonal Polynomials appear in the study of special solutions to the dispersive non- 
linear partial differential equation called the Degasperis-Procesi (DP) equation, as well as in certain 
two-matrix random matrix models. Another context in which such biorthogonal polynomials play a 
role is the cubic string; a third order ODE boundary value problem — /"' = zgf which is a generaliza- 
tion of the inhomogeneous string problem studied by M.G. Krein. A general class of such boundary 
value problems going beyond the original cubic string problem associated with the DP equation is 
discussed under the assumption that the source of inhomogeneity 51 is a discrete measure. It is shown 
that by a suitable choice of a generalized Fourier transform associated to these boundary value prob- 
lems one can establish a Parseval type identity which aligns Cauchy biorthogonal polynomials with 
certain natural orthogonal systems on Lg. 
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1 An ordinary inhomogeneous string, the Degasperis-Procesi 
equation and a cubic string 

1.1 DP equation 

The vibrating string is one of the most fundamental physical phenomena, whose mathematical foundations 
go back at least to the times of D. Bernoulli with important contributions of J.B. Fourier, B. Riemann, 
S. Sobolev and M.G. Krein and it appears that the heated discussion of more refined aspects of the 
problem has not subsided. Of multiple reasons for the special role of the vibrating string two stand 
out: 1) this is the simplest problem with infinitely many degrees of freedom and as such it inevitably 
leads to questions of analysis on function spaces, in particular the eigenvalue problems for infinitely large 
matrices; 2) the vibrating string problem is a prototype of hyperbolic problems and plays an important 
role in basic mathematical education. In its simplest, one dimensional, form, one studies the partial 
differential equation 

— Utt - U^xx =0, < X < 1, 

for the amplitude of the string u{x,t). The coefficient ^ is proportional to the mass density p of the 
string and in principle it does not have to be homogeneous (constant) in x. Thus, in general, = 
where T represents the tension in the string and p{x) is the linear mass density of the string. Then the 
fundamental modes, u{x,t) = v{x) coa{Lut), are given by solving: 

-Vxx = —p{x)v, < X < 1, 

subject to some boundary conditions, corresponding physically to how the string is tied. The most 
common are the Dirichlet boundary conditions, w(0) = v{l) ~ (the string is tied on both ends), and 
the Neumann conditions, Vx{0) = ^^(l) = (the end points can slide in the vertical direction only). 
Once equipped with the boundary conditions, say the Dirichlet conditions, the problem is turned into an 
eigenvalue problem of the Sturm-Liouville type: 

-Vxx^zp{x)v, v{0)=v{l)=0, 0<X<1 (1-1) 

where z = Such an inhomogeneous string problem was studied in an influential work by M.G. Krein 
[1] in the 1950s as a generalization of Stieltjes' theory of analytic continued fractions [2]. 
The present paper addresses several aspects of Krein's theory for a third order equation: 

-(jjxxx = Zp(j), < X < 1 

subject to some boundary conditions specified later. We will refer to this equation as the cubic string. 
Since any third order problem is non-self adjoint one can seriously doubt any applicability of this type of 
equation to physical systems. Below we briefiy sketch how the cubic string has arisen in a study of certain 
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nonlinear partial differential equation modeling weekly dispersive waves. The equation in question is the 
Degasperis-Procesi (DP) equation [3]: 

Ut - U^xt + 4:UUj: = SUxU^x + UU^xx, {x,t)eR^, (1-2) 

where u{x, t) is the wave hight at x and time t. The DP equation admits a Lax formulation, first proposed 
in [4] . which means that the DP equation follows from the compatibility condition for the the system 

{dx — <9^)V' = zmip, (l-3a) 
dti'^[z-\l^dl) + Ux-udx]iJ. (l-3b) 

where z G C. In general m can be a Radon measure, for example m can be a discrete measure m = 
Sr=i '^i^xi, and, as a result, the equation (1-3) will be assumed to hold in the sense of distributions. It 
is helpful to bring (l-3a) to its canonical form [5]. To this end one performs a Liouville transformation on 
(l-3a). This is fully explained in [5] and here we only state the essential results leading to the appearance 
of the cubic string boundary value problem. 



e^tanh-, V(x) = -J^, (1-4) 



Lemma 1.1. Under the change of variables 

I- '-(^'^^ 

the DP spectral problem (l-3a) is equivalent to the cubic string problem 

-</'«c(0= 2.9(0 0(0 /oree(-l,l), 

,/.(-l)=<^a-l)=0, (1-5) 
0(1) =0, 

where ^ 

5(0 = H^)- (1-6) 
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In the discrete case, when m{x) = 2 to^ (J^^ . , equation (1-6) should be interpreted as 

n „ 

5(0 = 2^5»^?i: where ^j=tanh— , gi = —2- (1-') 

.=1 ^ (1 - 0') 

To solve the original DP equation requires solving the spectral and inverse spectral problem as ex- 
plained in detail in [5] for the case of the finite discrete measure m. Our goal in the remainder of the 
paper is to present a few essential aspects of the spectral problem associated to the cubic string for 
a variety of boundary conditions with the view towards explaining the role of a novel class of Cauchy 
biorthogonal polynomials introduced in [6, 7]. These polynomials can be defined for arbitrary positive 
measures da, d(3 supported on the positive part of the real axis M+ provided that all (Cauchy) bimoments 
lij ~ J ^r^da{x)dP{y) are finite. One then defines polynomials qj{y),pj{x) of degree j = 0, 1, . . . satisfy- 
ing the biorthogonality condition J ^^J^y^^'^ da{x)df3{y) ~ Sjk- We refer to these as Cauchy biorthogonal 
polynomials. 
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2 Discrete cubic strings 



Wc slightly generalize the cubic string discussed in the previous section in connection with the DP equa- 
tion. We recall that in an ordinary string problem different boundary conditions correspond to different 
ways of tying down the ends of the string. For us, different boundary conditions will eventually lead to 
different spectral measures with respect to which we will define biorthogonal polynomials. However, in 
all cases discussed by us the spectrum is positive and simple, as one would expect from any vibrating 
system. 

Definition 2.1. The cubic string boundary value problems: 

-M0=^^9{0m, 0<e<l, 0(0) - 0^(0) = 0, 

(2-1) 

Type (DP case ): =0, Type 1 : (t)^{l) = Type 2 : (t>^^{l) ^ 
Remark 2.1. For simplicity we have adjusted the length of the string; it is now 1 rather than 2. 

We are only interested in the case where the mass distribution consists of a finite collection of point- 
masses: 

n 

giO = '^9^ . where < gi, < < 6 < ■ ■ ■ < 6i < 1- (2-2) 

i=l 

We will consider all three boundary value problems mentioned above with this mass distribution, as well 
as one degenerate case in which the last mass is placed at 1 (i.e. ^„ = 1: in the latter case we take the 
right hand limit to compute the derivatives of / at 1. Moreover, for that case, we consider only the BVP 
of type 2.) 

We will collectively refer to all these cases as the discrete cubic string problem. 
We will also use an accompanying initial value problem, 

0(0) =0^(0) =0,0^^(0)^1. (2-3) 

The boundary value problems in Definition 2.1 are not self-adjoint and the adjoint boundary value 
problems play an important role. 

Definition 2.2. The adjoint cubic string boundary value problems: 

reedo^^giorio, o<e<i, (2-4) 

Type 0: 0*(1) = 0|(1) = 0,^(0) = 0, Type 1: 0*(1) = 0|^(1) = O,0*(O) = 0, 
Type 2: 0|(1) = 0|^(1) = O,0*(O) = 0. 

The corresponding initial value problems are: 
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Definition 2.3. The adjoint cubic string initial value problems: 



< e < 1 



(2-5) 



with nonzero initial values: Type 0: 0|^(1) = 1, Type 1: <?I>|(1) = 1, Type 2: 0*(1) = 1- 

Since all three boundary value problems for satisfy the same initial value problem we will use one 
letter, namely 0, to denote the solution. However, we will attach an index a = 0, 1,2 to 0* to indicate 
the type of the EVP; for example, (jj^ wiU refer to the BVP/IVP of Type etc. 

In the process of integration by parts of expressions like /ccc(C)^(C)'^C one identifies the relevant 
bilinear symmetric form: 

Definition 2.4. Given any twice differentiable /, h the bilinear concomitant is defined as the bilinear 
form: 



This bilinear symmetric form induces a bilinear symmetric form (denoted also by B) on triples F'^ := 
ifJiJii), namely 



Furthermore, we also define a natural space associated with g, denoted i^[0, l]g, equipped with the 



we write {f{z)^ g{X))g to display this dependence. 

Since all initial value problems 2-3, 2.3 can be solved for arbitrary z S C, (/) and 0* are functions of 
the spectral parameter z. The following theorem establishes a relation between these two functions in 
terms of their concomitant and the relevant boundary value problem. The customary notation: /(^)|J = 
/(I) — /(O) is used in the statement below. 

Lemma 2.1. Suppose 0(^; z) and 0*(C; A) are solutions to the IVPs 2-3, 2.3 with spectral parameters z, 
A respectively. Then 

1. Type 0: the spectrum is determined by the zeros of (f){l, z) ~ and 



B{f, h){0 = feeh^fih + fhi. 



(2-6) 




(2-7) 



inner product: (/, h)g = fiO^iOdiO Whenever / or g depend on the spectral variables z and A 



sm-, ^),m. ml = m-, a) - 0(1; z)^{z- xmz),c^i{\))g. 



(2-8) 



2. Type 1: the spectrum is determined by the zeros of(f)^{l,z) ~ and 



sm-, ml = 0t (O; A) + (/.^(l; z) = {z~ A)(0(z), ^/.^ (A)),. 



(2-9) 
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3. Type 2: the spectrum is determined by the zeros 0/0^^(1, z) = and 

- ^),m. ml - </'2(o; a) - ^^(i; ^) = - a) w^), ^sIa))^. (2-10) 

In addition, 

- Bm- z),(b{C, ml = (2 + A)(0(z), 0(A))g. (2-11) 
Proof. Indeed (2-3) and two integrations by parts imply that 

^)m: A) = - i?(</', ^*a)\l+ f zWa,i^^{e, A) = 

Jo 

z C m^Wa{i\m{i)di. 

Jo 

Consequently, using equation (2-5) we obtain: 

-B{cp, 0:)|i = {z-\) f (/.(C; z)</>:(e; A)g(0 d^, 



^0 

which in view of the initial conditions implies the claim. A similar computation works for the second 
identity. □ 

By specializing A = z in the lemma above one readily obtains: 

Corollary 2.1. (j) and 0* satisfy the following relations: 

1. Case 0: 0^(0; z) = 0(1; z). Case 1: 0^(0; z) ^ -0^(1; z), Case 2: 4>*2(0]z) = ^^^(l;^)- 

2. Case 0: -0.(1; z) = (0(z), ^^(z))^. Case 1: 0^.(1; z) = (0(z), 0t(z))g. Case 2: -0c5.(l;z) = 
(0(z),0*(z)),. 

We give below a complete characterization of the spectra and the corresponding eigenfunctions for all 
three BVPs. We also select certain combinations of EVPs which reveal the origin of the relevance of the 
Cauchy kernel ^ji^ to the spectral theory of the cubic string. 

Theorem 2.1. Consider a cubic string with a finite measure g as in (2-2). 

1. Let Zaj denote the eigenvalues of the EVP of Type a ~ 0,1,2. In each of the three cases, the 
spectrum is positive and simple. 

2. For any pair of BVPs of Type 0,1,2 the spectra are interlaced in the following order: 

Z2,j < zi,j < zqj, j = l,...n 

3. For each a = 0, 1,2 the eigenfunctions Zaj) = 0oj(C) ^'"6 linearly independent. 



4-. For the following combinations of BVPs, {4>a,i,4'b,j)g factorizes: 
(a) Type 00: 



(0o.,<^o.-).- ^°-^^^^-'^^^^ (2-12) 



(h) Type 01 

{(f>0,i,<l)l.,j)g 



0o,^$g(l)0ij(l-) 
zo,i + zi.j 



(c) Type 12 



(d) Type 22 





where 




and 



(2-13) 



iM,<A2.). = -^i^^^^^^^^ (2-14) 



Zl,i + Z2J 



h,), = (2_i5) 

Z2A + Z2.j 



Proof. It is easy to check (see Section 4.1 in [5]) that 

= L„Gn{z) L„_i G„_i(z) • • • Gi(z) io I I . (2-16) 



lk=^k+i-^k, 6 = 0, Cn+i^l. (2-17) 



We prove the statement about the spectra by using the results obtained in [5]. By Theorem 3.5 in [5] 
0(1; z) has n distinct positive zeros and so do 0^(1; z) and (^^^(1; z) (denoted there (l)y, 4>yy). The second 
statement foUows if one observes that and '^^''(^i'^^ are strictly positive on the spectrum of Type 

being the residues of ^j:^, t^^i^^)'' respectively (Theorem 3.5 and Theorem 3.15 in [5]). It follows then 
that both 0^(1; z) and 0jj(l;z) change signs n times, hence all three spectra are simple. Furthermore, 
0^(1; z) > 0, (j)^^{l;z) > for z < 0, so the zeros of z) and (j)^^{l;z) are strictly positive. Since 
0^(1; z) and (jj^^i^'y z) change sign on every interval between two consecutive zeros of 0(1; z) the spectra 
of Type and 1, as well as and 2, interlace. To see that the spectrum of Type 1 interlaces with the 
spectrum of Type 2 we proceed as follows. By (2-11), after evaluating at z = Z2,i, A — ~Z2,i, we obtain 

0e(i; ^2.j)0c(i; ~z2,i) = 055(1; -Z2.J0(1; Z2,,), 
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which gives 

05(l;z2,) = MiiZf^0(l;,,,^) (2-18) 

For z > 0, sgn(0j(l; — z)) ~ sgn((/)^^(l; — z)) = +1 because both are strictly positive there. Since the 
zeros of 0(1; z) interlace with the zeros of ^^^(l; z), sgn((/)(l; Z2a)) alternates, which in turn implies that 
sgn((/)^(l; Z2,i)) alternates as well. Thus the zeros of 0^(1; z) interlace with the zeros of 0^^(l;z). 

The relative position of the spectra of the three types is best inferred from the fact that on the first 
eigenvalue zo,i of Type 0, 0^(1; zq^) and ^^^(l; zo,i) are both negative since and '^^^l^i'^^ are strictly 

positive on the spectrum of Type 0. Thus the spectra of Type 1 and 2 are shifted to the left relative to 
the spectrum of Type 0. In particular, 0(1, Z2,i) > and so is 0^(1, Z2.1) by (2-18). Thus, at least the first 
zero of 0JJ occurs to the left of the zeros of 0^ and 0. So Z2,i < zi i < zo,i. Suppose this holds for the 
{j — l)st eigenvalues. Then we know that both zo,j_i < zi j < zqj and zoj_i < Z2j- < zqj. If zi j < Z2j 
then Z2j_i < < zij < Z2j , thus contradicting that the spectra of type 1 and 2 interlace. 

As for the linear independence we observe that the cubic string boundary value problem 2.1 can be 
written as an integral equation: 

0(^; z) = z [ G(^, t)0(t; z)g(T) dr, (2-19) 

JQ 

where G(^,r) is the Green function satisfying the boundary conditions of 2.1. Then the linear indepen- 
dence of eigenfunctions corresponding to distinct eigenvalues follows from the general result about the 
eigenfunctions of a linear operator. 

Finally, the statements about the eigenfunctions follow immediately from equation (2-11) after setting 

Z — Za^ii A — ZbJ' 

□ 

We now briefly analyze the degenerate case with the mass m„ at the end point a;„ = 1. We do it only 
to illustrate that even though the spectrum degenerates in this case the overall conclusions hold. 

Theorem 2.2. Let us consider a cubic string with a finite measure g as in (2-2) with x„ = 1. Then 

1. the spectra of Type 0, 1, and 2 are positive, simple and their cardinalities are n — 1 for Type and 
Type 1 and n for Type 2, 

2. the spectra interlace 

< Z2,l < Z14 < Zo,l < • • ■ < Z2,„-l < 2l,„-l < Zo,„-l < Z2,„, 

3. the eigenfunctions 0(^;z2j) :~ 02.j(C) satisfy 

/ , , N 02,»,^(1)02,m(1) 

(02,j,02j)3 = ■ , (2-20) 

Z2d + Z2J 

and they are linearly independent. 
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Remark 2.2. We notice that deg(t){l; z) = deg0^(l; z) = n—1 while dcg0^j(l; z) ~ n. This is in contrast 
to the previous cases with all positions Xi, . . . ,Xn inside the interval [0, 1] for which all polynomials have 
the same degree n. 

Proof. The spectrum of Type 2 is clearly given by the zeros of <f>^^{l;z). Let us first consider the case 
when m„ is placed slightly to the left of the point 1. Thus, initially, /„ > (see (2-16)). By Theorem 2.1 
0(1; z) has n distinct positive zeros and so does and they interlace. Wc subsequently take the limit 
In in the above formulas. We will use the same letters for the limits to ease the notation. By simple 
perturbation argument, z„^o ^ cxd- Since z = is not in the spectrum, Z2.1 has to stay away from 0. This 
shows that the spectrum is positive. Furthermore, in the limit zo,ij • ■ • j ^o.n-i approach simple zeros of 
the BVP of type for n — 1 masses. Indeed, using (2-16) with In = there, we obtain: 

0(1 + 0;z)^ 0(1 - 0; z), 0^(1 + 0; z) = 0^(1 - 0; z), 0^^(1 + 0; z) = -zm„0(l - 0; z) + 05^(1 - 0; z), 

(2-21) 

where 1 ± refers to the right hand, the left hand limit at 1 respectively. To see that the spectrum 
is simple we observe that if in the limit two successive eigenvalues coalesce, namely Z2,i = Z2^i+i, then 
necessarily Z2.i ~ zqj because of the interlacing property. However now equation (2-21) implies that 
0(1 — 0; ZQ i) = 0^5(1 — 0; ZQ i) ~ which contradicts Theorem 2.1 for the BVP of Type for n — 1 masses. 
Thus, the zeros of 0(1, z) and 0^^(1 + 0, z) interlace and we have 

< Z2,l < Z04 < ■ • • < Z2,„-l < Zo,„-l < Z2,„. 

Likewise, for the spectrum of Type 1, Zn^i — > 00 and the remaining roots interlace according to the 
pattern valid for n—1 masses. 

To prove the statement about the eigenfunctions we use (2-11) and after setting z = z^, A = zj in that 
formula we obtain the required identity. The linear independence is proven by the same type of argument 
as in the proof of Theorem 2.1. □ 

We immediately have several results about the adjoint cubic string 2.2. 

Corollary 2.2. Given a discrete finite measure g 

1. for each of the three types of the BVPs the adjoint cubic string (Definition 2.2) and the cubic string 
(Definition 2.1) have identical spectra. 

2. the families of eigenfunctions {0a. j} and {0aj} are biorthogonal, that is: 

(0a,i, 0*j)g = whenever i^j- (2-22) 
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3. For i = i, 

r-0,(l;zo,,:) 7^0, a = 

= <^ 0^,(1; zi,,) ^0, a = l (2-23) 

[-045,(1; zz.z) 7^0, a=2 

holds. 

Proof. The first equality in Corollary 2.1 implies that the spectra of the cubic string and its adjoint are 
identical. The biorthogonality follows immediately from equations (2-8), (2-9) and (2-10). For i = j we 
use Corollary 2.1. Finally, since the spectrum is simple the required derivatives with respect to z are 
nonzero. □ 

Wc conclude this section with the definition and some fundamental properties of two important 
functions which play a significant role in the theory (Section 6 in [6]). 

Definition 2.5. The following functions are called Weyl functions for their respective BVPs: 
Type : W{z) := Z{z) = Type 1 : W{z) := Z[z) = 

Type2 ■.W{z):=^P^^, Z{z) - ^^''^ 



055(1; z)' 055(1; z)' 

Remark 2.3. The definition of the Weyl functions for the EVP of type 2 in the degenerate case is 
identical to the one given above for the B VP of type 2. 

The Weyl functions W and Z arc not independent, they arc related by an identity which was originally 
formulated for the DP peakons in [5] . As an example wc formulate such an identity for the BVP of type 
2 (both the degenerate as well as the nondegenerate case). 

Lemma 2.2. Consider the BVP of type 2. Then the corresponding Weyl functions satisfy: 

W{z)W{~z) + Z{z) + Z{-z) = Q (2-24) 

Proof. By formula (2-11) 

B(0(e;z),0(e;-z))|i=O, (2-25) 
which, when written out explicitly, gives the identity: 

055(1; z)0(l; -z) - 05(1; z)05(l; -z) + 055(1; -z)0(l; z) = 0. (2-26) 

Upon dividing the last equation by 055(1; z)055(l; — z) we obtain the claim. □ 

We state now the spectral representation theorem for W{z) and Z{z) for the BVP of Type 2 in 
the degenerate case, the remaining cases being merely variations of this, most transparent case. One 
recognizes again the presence of the Cauchy kernel in the spectral representation of Z{z). 
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Theorem 2.3. Consider the BVP of type 2 (degenerate case). Then the Weyl functions W and Z have 
the following (Stieltjes) integral representations: 

W{z) = / -^d(3{y). Z{z) - / i-^d/3(2/)d/3(x), (2-27) 

J z-y J {z-y){x + y) 

where d(3 = YJi^i h5z2,^ and bi = ^^'^^(l'.^^''.) > 0. 

Proof. Since z), ^^^(l; z) have simple, interlacing zeros, and deg z) = n — l while deg (/'{{(I; z) ~ 
n, W{z) admits a partial fraction decomposition with simple factors: 

br 



W{z)^Y. 



. Z~ Z2., 
I— I 



where, by the residue calculus, bi = t^^tt^-^- Moreover the 6iS are all of the same sign because the 



zeros of 0^(1; z) and ^^^(l; z) interlace and, consequently, it suffices to check the sign of 0^*^^ at the 

first zero Z2.\- By Theorem 2.1 sgn((/)^ (1; Z2,i) = 1, and thus &i > because on the first zero (^i^^f^z must 
be negative. Consequently, all bi > 0. 

Likewise, Z{z) admits a partial fraction decomposition: 



Ci 



- Z - Zl.i 

and again, it follows from the second item in Theorem 2.1 that q > 0. Finally, by residue calculus, it 
follows from Lemma 2.2 that 

Ebibj 
, 

which proves the integral representation for Z(z). 

□ 

3 Generalized Fourier transform and biorthogonality 

Since (/)a,i are linearly independent we can decompose any / G ^^[0, 1] in the basis of {0a, i} and use the 
dual family {(^* j} to compute the coefficients in the expansion: 

For each pair a,b for which {(t>a,i,4'b,j)g factorizes (item 3 in Theorem 2.1) we define two (finite dimen- 
sional) Hilbert spaces Ha := i^(]R, da) and Hp := L^(R, d/?) where the measures da and d(3 are chosen 
as follows: 1) using Theorem 2.1 split the numerator of {4'a.ii 4'b.j)g , 2) perform the partial fraction 
decomposition of the Weyl fimctions with numerators matching the factors in {4'a,i, 4'b.j)g , 3) represent 
the partial fraction decompositions as Stieltjes' transforms of the respective measures. 
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Example 3.1. For Type 00, item 3 in Theorem 2.1 states: {(po.i, (j>Qj)g = ■ Hence the 

numerator splits into (/)o,,;.^(l) and (j3oj.^{l) and the corresponding Weyl function will be W{z) = '^^li-z) 
(taken twice) with partial fraction decomposition W{z) = J j^da{x) where da = dP = ^ "(I'-zo ) ^^<> ' 
For Type 01, item 3 in Theorem 2.1 states {<po.iT 4'i.j)g = — and the numerator splits into 

0o.i,{{(l) and 01 J (1) which, in turn, match with Z{z) for Type 1 and W{z) for Type 2 in the Definition 
2.5 . The resulting partial fractions decompositions Z{z) = j j^da{x), W{z) ~ I T^d.p{y) imply 

In summary, for every pair {a, b} of BVPs for the cubic string we associate two Hilbert spaces Ha 
and Hp with the pairing Ha x Hp — + C: 



Definition 3.1. 



{p\q) = f ^^^^da{xmy), p € Ha, g G Hp. 
J x + y 



We now introduce a family of generahzed Fourier transforms adapted to each of the three types of 
BVPs 

Definition 3.2. Given f e Ll[0, 1] and a = 0, 1, 2 



fa{z):^{-lT m;^)f{i)9{m- (3-1) 
Jo 

Remark 3.1. Observe that fa{z) = ((— 1)°0* (z), and, in particular, fa{za.i) = ii—^)'''4>a.ijf)g 
whenever z equals to one of the points of the spectrum of (2-1). 

Remark 3.2. A map of this type was introduced in the context of the inhomogeneous string problem by 

LS. Kac and M.G. Krein in [8] as a generalization of the Fourier transform. A similar map known as a 
distorted Fourier transform is commonly used in Quantum Mechanics. 

The main property of this map is captured in the foUowing theorem. 

Tiieorem 3.1. For every pair a, b for which {(f>a.i, 4'b,j)g factorizes (item 3 in Theorem 2.1) the generalized 
Fourier transforms satisfy FarsevaVs identity, that is, for every f,h^ ^V^^^ ^'^ 

U,h)g^{fa\h) = {ha\h). (3-2) 

Proof. Let us fix a,b for which {(j)a,i,4>b,j)g factorizes. Consider two functions f,h ^ Lg[0, 1]. Writing 
their expansions in the bases {0a, i}, {4'b,i} respectively, we obtain 

. i(t>hJ)g , , (<^bj'Ms , 
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Hence their inner product reads: 

Applying now item 3 from Theorem 2.1 as weU as item 2 from Lemma 2.1 we obtain 

TypeOO :(/, h), = J] I^.T '\ ^^^^^ ^"'^^^^^^ ^"^^ 

^ 02(1; zo,i) 0^(1; zoj) ^o,* + zo,j 

^ 02(1; Zo,i) 0jz(l;2:ij) ^O.i + 



'?!'Cz(i;zi,j) '/'ccz(i;2:2j) + 22, 
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^ '?'C5z(l; 22,i) <P??2(1; 22j) Z2,i + 22j 

We now define the weights hj^a^ generating the measures d/3 = Bjd^j^ da ~ ^j^zj respectively, 
as residues of W^s or Zs: 



<Pz{1,zoa) 05^(1; zij) 



j(l:2l,»)' ' 055^(1, Z2j)' 

Type 22 : A, = ^^fe^, = -^^^^^'''^^ 



and thus indeed 



(l)iizil,Z2,i)' ^ 055^(1, 22,^)' 



- ■ Z(2 i -f- 2:^ 



Thus, in view of Remark 3.1, we obtain 



□ 



Remark 3.3. Expanding an arbitrary f £ L'^[0, 1] 

„ _ ('^a,i' f)g , 

allows one to conclude that 

plays the role of the Dirac delta on Lg[0, 1]. 
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It is now elementary to find the inverse Fourier transforms 

Lemma 3.1. Consider the BVP of type a. Let {za.i\ he the corresponding spectrum and let dva — X)i ^Za. i 
he an associated measure. Then the inverse generalized Fourier transform of type a is given hy 

'-'>i "-'WW,'''''''' <'-^' 

Proof. This is a direct computation: 



, [4>ni^Va,i)g Jo 



□ 



There are in general two measures associated with each type of the BVP, one generated by W the 
other by Z. One can use either one of the them instead of the measure di^. We give as an example the 
relevant statement for the the case of the BVP of Type 2, both the degenerate and the nondegenerate 
one. 

Lemma 3.2. The inverse generalized Fourier transform of Type 2 is given hy 
Proof. From the definition of d/3 given in Theorem 2.3 we see that 

[ f ( \ ^^^''^) ^ni \ \ 'l^^-AO , f f \ 'f'^AO / <?!'j(i;^2,*) s 

J M^'^^) i ?!'?(l;^2,0 ^ 05(1; Z2,,;) (/)55^(1,Z2^,) 

^ 0«(l;22,j) (l)iiz{l\Z2,i) ^ (02,j>2,j)ff 

where in the last two steps we used Remark 3.1 and equation (2-23) respectively. Thus 
J M^'^^) Jo i [<f>h^(l'2Ag Jo 

□ 

The generalized Fourier transform can be used to give the following interpretation of the Cauchy 
biorthogonal polynomials associated to discrete, finite positive measures da, dp. 
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Theorem 3.2. Biorthogonal polynomials Let us consider the sequence Xj '■— X(5„_3-£,{„_j+e) of 
indicator functions enclosing points £,n-j of the support of the measure g appearing in the cubic string 
equation (2-1), with e small enough to ensure non overlapping supports. Furthermore, let us consider the 
generalized Fourier transforms for Types a,b as in Theorem 3.1 and define: 

Pj{x) ■= Xajix), qjiy) := X&j(y). 
Then, {pj, qk} are (un-normalized) Cauchy biorthogonal polynomials, that is < Xa.i\Xb,j >= 0, i j and 

Proof. Both Xa,iix) and Xbjiu) are polynomials in x, y respectively, whose degrees are deg Xa.j{x) = 
degXb,j{y) = j by (2-5). The biorthogonality follows from Theorem 3.1. □ 
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